Long-Range Dynamic Correlations in Confined Suspensions 
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Hydrodynamic interactions between particles confined in a liquid-filled linear channel are known 
to be screened beyond a distance comparable to the channel width. Using a simple analytical 
theory and lattice-Boltzmann simulations, we show that the hydrodynamic screening is qualitatively 
modified when the time-dependent response and finite compressibility of the host liquid are taken 
into account. Diffusive compression modes in the confined liquid cause the particles to have velocity 
correlations of unbounded range, whose amplitude decays with time only as t~ 3 ^ 2 . 
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As particles move through a liquid, the flows that they 
create correlate their motions. These hydrodynamic in- 
teractions play a central role in the dynamics of partic- 
ulate liquids In unbounded liquids the correlation 
is dominantly mediated by diffusive shear modes (i.e., 
diffusion of transverse momentum) as accounted for by 
the Navier- Stokes equation for an incompressible liquid. 
Consequently, the perturbation caused by the forced mo- 
tion of a particle in the liquid can be represented at large 
distances as emanating from a momentum monopolc. 
The resulting flow velocity field decays with distance as 
1/r, leading to long-range hydrodynamic interaction be- 
tween particles and the well studied phenomenon of long 
time tails in their velocity correlation functions 0] ■ 

When the suspension is confined by rigid boundaries, 
the hydrodynamic interaction is qualitatively different 
Q . Because of momentum transport into the boundaries 
the liquid no longer conserves momentum, and the ampli- 
tude of shear modes, developing in the liquid in response 
to particle motion, decays exponentially with distance. 
The correlations are now dominantly mediated by com- 
pression modes, and the perturbation due to particle mo- 
tion can be represented at large distances as arising from 
an effective mass dipole. When the suspension is confined 
between two plates (quasi-two-dimensional suspension), 
the resulting flow field remains long-ranged, decaying as 
1/r 2 (as for a mass dipole in two dimensions) QJ. How- 
ever, since in one dimension (ID) a steady dipole has 
a vanishing effect, the steady flow due to local particle 
motion in a linear channel has no far field and decays 
exponentially with r/h, h being the channel width 
This hydrodynamic screening in pores and channels is 
well documented (e.g., 0, 0) and widely employed (for 
instance, in studying the dynamics of confined polymers 

!)• 

The analysis above assumes that the host liquid is at 
steady state and is incompressible, thus omitting lon- 
gitudinal sound modes. Since the Reynolds number of 
relevant colloidal systems is typically smaller than 10 ~ 6 , 
and the sound velocity in the liquid, c s , is of order 10 
m/s, both assumptions seem safely valid. Yet, as was 
first pointed out in Rcf. 9] based on simulations, and 



further established analytically lfl U|, confinement by 
rigid boundaries strongly affects the sound modes in a 
compressible liquid. These modes change from under- 
damped propagation with velocity c s to overdamped dif- 
fusion with diffusivity D s ~ c 2 /i 2 /^, where v is the kine- 
matic viscosity of the liquid. The resulting sound dif- 
fusion is fast, D s ~ 1 m 2 /s for water in micron-scale 
confinement, compared to v ~ 10~ 6 m 2 /s for the diffu- 
sion coefficient of shear modes, yet the latter are sup- 
pressed by the boundaries, leaving the diffusive compres- 
sion modes as the sole mechanism for hydrodynamic in- 
teraction at interparticle distances r > h. Furthermore, 
similar to the appearance of long time tails due to shear- 
stress diffusion in unbounded suspensions (2j, the diffu- 
sive nature of compression modes in confined liquids leads 
to a (negative) long time tail in the velocity autocorrela- 
tion function of a suspended particle [9j. In this Letter 
we show that those diffusive compression modes have far 
reaching implications for the dynamic correlations be- 
tween distant particles in a channel. 

We employ a simple phenomenological approach, in 
which only the unconfined dimensions (in the current 
case the single dimension along the channel) are retained, 
whereas all the details in the confined dimensions, such as 
the cross-section of the channel and the precise boundary 
conditions at its walls, are averaged over and reduced to 
an effective friction term. Such a lubrication approach 
[lH leads to tremendous simplification of the analysis 
while yielding the correct physics at distances larger than 
the confinement width. It was successfully applied in the 
pas t to various problems of confined hydrodynamics (e.g., 
[131 ] ) and, in particular, was shown to correctly reproduce 
the diffusive sound modes in confinement [9I lll|. 

We support the calculations by three-dimensional 
lattice-Boltzmann simulations [Hj]. In the simulations 
spherical particles of diameter a move along the axis of a 
channel of length L and a square cross-section of side h. 
The channel is filled with an ambient fluid of mass density 
po, shear viscosity 77, bulk viscosity £, and sound veloc- 
ity c s . The channel and particles are rigid, and a no-slip 
boundary condition is imposed at their surfaces using the 
second-order bounce-back scheme [3]. Periodic bound- 
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ary conditions are imposed at the edges of the channel, 
though L is taken sufficiently large to render finite-length 
effects negligible during the simulated time. In terms of 
lattice spacings and simulation time steps the parameter 
values are h = 15, a = 10, po — 1, i] — 1/6, ( — 1/9, and 
c s = l/\/3- To present the simulation results in a way 
relevant to real systems we shall use as units of length, 
time, and diffusivity, respectively, the channel width h, 
the time it takes shear modes to diffuse to the bound- 
aries, t = h 2 jv {y — i]/po), and the self-diffusion coeffi- 
cient of an unbounded particle, Do = k B T '/ '{iir-qa) , k B T 
being the thermal energy. For particles of a = 1 pm in a 
water-filled channel of h — 1.5 /an at room temperature, 
we have r ~ 2 ps and D ~ 0.4 pm 2 /s. 

Within the simplified theory the linearized ID hydro- 
dynamic equations, corresponding to a channel, read 



p ii = -p' + (477/3 + C)u" - (a^/h^u + ((3/h 2 )f, 



p = -p u', p = c 2 p. 



(1) 



Being one-dimensional, Eq. (fTJ) describes solely longitu- 
dinal flow. In it a dot denotes a derivative with respect 
to time t, and a prime is a derivative with respect to dis- 
tance x along the channel. The fields p(x,t), p(x,t), and 
u(x,t) are, respectively, the perturbations in liquid mass 
density, pressure, and velocity about po, Po, and 0, in re- 
sponse to the force density (per unit length) f(x,t). We 
have introduced two dimensionless coefficients — a, de- 
termining the strength of the effective friction caused by 
the boundaries, and (3, equal to the ratio between h 2 and 
the cross-sectional area of the channel. For the simulated 
square channel j3 = 1, and we find a ~ 28 by a fitting 
procedure independent of the following analysis fl5j . 

The velocity Green's function, u(x,t) = G(x,t), 
obtained from Eq. (JlJ upon substituting f(x,t) = 
5(x)5(t), is readily calculated in Fourier space [G(q,u) — 
JJ dxdte-^-^G&t)} as 
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(r) + ic 2 p /uj)q 2 + arj/h 2 — ipoUJ ' 
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where rj = 4^/3+C. Inverting the spatial coordinate back 
to real space, we obtain 

G(x,u) = A(w)e- |x|/A(w) 

A(u) = [(3/(2h 2 )} [(fj + ic 2 s p /u)(a V /h 2 - ip uj)] ^ 

^ ^ ^arj/h 2 — ipoUJ^) 
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The function G(x, t) gives the flow velocity in the channel 
as a function of position and time in response to a point 
impulse applied at the origin at t = 0. If two particles 
are positioned inside the channel at sufficiently large mu- 
tual distance, x ^> a, the same function, up to a dimen- 
sionless factor, -f(a/h), gives the velocity of one particle 
in response to a unit momentum imparted to the other. 



Hence, the velocity cross-correlation function of the par- 
ticle pair is given, according to the fluctuation-dissipation 
theorem, by C{x,t) = (Fi(0)F 2 (t))(a;) = jk B TG(x,t) 
The coupling diffusion coefficient, characterizing 
the correlated Brownian motion of the pair, is then ob- 
tained from the appropriate Green-Kubo relation 0, 
D c (x) = J dtC(x,t) = 7feTG(j;,a) = 0). 

Two important observations readily follow from Eq. 
([3|). First, since A{u = 0) = 0, the coupling diffusion 
coefficient vanishes, D c (x) = 0. This is the manifesta- 
tion of steady-state hydrodynamic screening in the chan- 
nel. The actual exponential decay of D c (x) arises from 
the screened transverse modes, which have been elimi- 
nated in the simplified ID calculation. It is clearly seen 
in the simulation results shown in Fig. [TJ Note that peri- 
odic boundary conditions have not been imposed in Eq. 
([3]), allowing a pressure difference to form between the 
edges of the channel 0, [ijj]. Second, as the diffusive re- 
sponse broadens, the correlation length increases indefi- 
nitely, A ~ cj -1 / 2 for small u>, while the correlation am- 
plitude decreases as A ~ a// 2 . This leads at long times 
to a correlation of unbounded range and slowly decaying 
negative amplitude, 



C(x,t^ oo) /(jk B T) 
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where D s = c 2 h /{av). To reach this distance- 
independent behavior the diffusive front must have al- 
ready passed the distance x, i.e., t > x 2 /D s . For x ~ 1 
cm in a micron-wide channel filled with water this re- 
quires t > 10~ 4 s, i.e., Eq. (0| holds practically at all 
times. (For liquids of higher viscosity, however, the 
distance-independent regime will take a longer time to 
establish.) The negative i~ 3 / 2 tail is confirmed by the 
simulation (Fig. [2j\ inset), from which we get for the 
square channel j(u/h = 2/3) ~ 2.61 17 1. 

In slightly more detail, the temporal correlation be- 
tween the velocities of two particles consists of three parts 
(Fig. [2]) . Because of the finite speed of sound there is a 
short incipient period of no correlation, which becomes 
longer with increasing interparticle distance. It is fol- 
lowed by a short peak of positive correlation upon the 
arrival of the signal and a subsequent long negative tail, 
which is the result of sound diffusion through the channel. 
As can be seen in Fig.[5j3, the simple theory qualitatively 
captures all three features. As the interparticle distance 
becomes larger than h, the quantitative agreement be- 
tween theory and simulation progressively improves. 

At steady state (w = 0) the correlation becomes ex- 
ponentially small in x/h (Fig. [T|). This is because the 
long-ranged temporal correlation arises from longitudi- 
nal modes, and these cannot produce a net displacement 
for any finite wavelength. (If they could, it would imply 
a nonuniform steady density of the liquid.) Thus, the ar- 
eas under the theoretical curves of Fig. 03 all vanish, and 
the remaining steady-state coupling is produced solely by 
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FIG. 1. Coupling diffusion coefficient of a particle pair in a 
channel as a function of interparticle distance, obtained from 
lattice-Boltzmann simulations. The coefficient is measured 
from the net displacement (time-integrated velocity) of one 
particle in response to a unit-momentum impulse imparted 
to the other. The particle diameter is a = h/1.5. 



the screened shear modes (the finite total areas under the 
simulated curves of Fig. [2j\). 

So far the way to make the unscreened longitudinal 
modes effective has been to include liquid compressibil- 
ity. Another specific but noteworthy example of an un- 
screened response is due to the uniform (q = 0) mode, 
which is unaffected by the incompressibility constraint. 
Indeed, from Eq. ((2J one has G(q ^ 0,u) = 0) = but 
G(q = 0, u) = 0) = p/{arj) > 0. Consider an infinite ar- 
ray of particles along the channel, separated by a distance 
d from each other, and examine the collective diffusion 
coefficient of the array, -D co i, corresponding to collective 
motion of all the particles. For d — > oo D co \ = D se \{, 
the self-diffusion coefficient of an isolated particle in the 
channel. Thus, D co \ = D co i-D sc i f = 2J2^=i D c{x = nd) 
characterizes the contribution to the collective motion 
from flow-induced correlations. With only transverse 
modes in mind, since each of the pair couplings D c (x) 
is screened (Fig. [1]), we should expect D co \ to decay 
exponentially with d/h. This behavior, expected from 
conventional hydrodynamic screening, is depicted in the 
inset of Fig. [31 However, considering the longitudinal 
flow of Eq. ([3]), we get 

oo 

D co i/(jk B T) = 2j2G(x = nd,Lo)\^ = (3/(a V d), (5) 
n=l 

i.e., a much slower decay of D co \ ~ 1/d. The appear- 
ance of 1/d in -D C oi, in place of the 1/a in Dq, implies 
that the particles, no matter how far apart, carry the 
liquid column in between them as they move collectively 
through the channel. This again reflects an interparticle 
coupling of unbounded range. As suggested by the ab- 
sence of po and c s in Eq. ([5]), this result is also obtainable 
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FIG. 2. (color online). Velocity cross-correlation function 
of a particle pair in a channel as a function of time. Dif- 
ferent curves correspond to different interparticle distances 
(from left to right): x/h = 1.33 (black), 3.33 (red), and 6.67 
(blue). The particle diameter is a = h/1.5. (A) Lattice- 
Boltzmann simulation results. The correlation is measured 
from the time-dependent velocity of one particle in response 
to a unit-momentum impulse imparted to the other. The inset 
shows the long-time behavior on a log- log scale, the dashed 
line indicating a slope of —3/2. (B) Results from the analyt- 
ical approximation [inversion of Eq. J5}]. 



from a simpler incompressible theory. In the q — case 
the incompressibility constraint, as well as the periodic 
boundary conditions, are trivially obeyed. There is no 
pressure difference between the edges of the channel and, 
consequently, each forced particle creates a uniform flow 
velocity ~ (rjL) -1 [f6j], which affects all others. Thus, the 
collective mobility of N particles is ~ NljjL)^ 1 — (ijd) -1 , 
in accord with Eq. (JSJ) - It follows that the same result 
is valid for a finite ring-like channel and for a disordered 
train of particles, whereby the factor 1/d in Eq. §5§ is 
replaced with the mean linear density of particles. 

The collective mode of an infinite particle array is eas- 
ily simulated by considering a single particle at the center 
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FIG. 3. Collective diffusion coefficient of an array of particles 
as a function of interparticle separation. Lattice-Boltzmann 
simulation results (solid) are presented alongside the analyt- 
ical result [Eq. ([5|; dashed]. The inset shows the collective 
coefficient as expected from the hydrodynamic-screening de- 
scription. The particle diameter is a = h/1.5. 



of a channel of length d and imposing periodic boundary 
conditions in the x direction. The collective mobility is 
then obtained from the net displacement (time-integrated 
velocity) of the particle in response to an initial impulse. 
As shown in Fig. [3j the simulation results are in good 
agreement with Eq. ([5]), exhibiting a slow 1/d decay. (No 
fitting parameters are used in Fig. [3] in addition to the 
already extracted values of a ~ 28 and 7 ~ 2.61.) 

In summary, the common view, according to which 
particles move through liquid-filled channels and pores 
in an essentially uncorrelated way, is quite misleading. 
Before steady state is reached, the motions of confined 
particles are in fact correlated over large distances and 
long times. The physical mechanism behind the correla- 
tions is fundamentally different from that in the uncon- 
fined case, as it originates from the confinement-affected 
longitudinal (rather than transverse) liquid response fljj . 
Consequently, the long-ranged correlations are negative 
(rather than positive) and do not give rise to correspond- 
ingly long-ranged steady correlations, so that at steady 
state hydrodynamic screening is recovered. These long- 
ranged temporal correlations may have important impli- 
cations, for example, for suspended particles in porous 
media or the dynamics of polymer translocation through 
narrow pores. Despite the small factor D s ~ c^ 1 
appearing in Eq. they should be observable when 
probing sufficiently short times. For instance, over a mil- 
lisecond (kHz frequency) Eq. (j4j predicts for particles in 
a micron-wide channel filled with water a velocity cross- 
correlation of order — 1CU 2 (/im/s) 2 , independent of par- 
ticle separation. In addition, the predicted strong effect 
that arises from the uniform longitudinal mode [Eq. ([5])] 
can be readily checked by measuring the collective mo- 



bility or diffusivity of a dilute particle assembly confined 
in a closed ring-like microfluidic channel. 
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In the simulation the velocity is measured on the chan- 
nel axis, whereas in the ID theory it represents a cross- 
sectional average; this affects the fitted value of 7. 
The exponent —3/2 in the velocity correlation function 
being identical to that in an unconfined liquid is a sheer 
coincidence. In the unconfined case it is equal to — D/2 
with dimensionality D — 3, whereas in the channel it is 
equal to -(D/2 + 1) with D = 1. 



